Entangling photons through nonlinear response of quantum wells to ultrashort pulses 
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We show that many-body correlations among excitons originating from the Pauli exclusion prin- 
ciple in a quantum well embedded inside a microcavity provide a possibility to produce pairs of 
entangled photons by ultrashort laser pulses with a yield of ~ 10~ 2 . The quantum-field theoretical 
two-particle density matrix in second quantization is used to calculate entanglement for arbitrary 
emission angles. Largest response can be expected at symmetric emission angles for resonances 
with the heavy-heavy and light-light two-exciton states with remarkably nontrivial dependence of 
entanglement on the emission angles and on the ellipticity parameters of the excitation. We show 
that the angle dependence can be tailored by means of the microcavity. Interestingly, the emitted 
entangled 2-photon states are always in a triplet state. 

o 

Current entangled-photon sources are mainly based on the parametric down-conversion (PDC) inside nonlinear 
crystals [Hli|, 

such as BBO crystals^, These sources suffer from two serious limitations. First, since the far off- 
resonant three-photon scattering contains two far off-resonant virtual states, the entangled-photon production yield 
is very iow[i, Gf, which limits the brightness of entangled-photon sources based on nonlinear crystals, leading to low 
signal-to-noise ratios and long measurement times [2]]. Second, PDC produces entangled photons with the twice as 
long wavelength as the pump photons, which limits the operating wavelength Therefore, the problem of alternative 



sources of entangled photons is of the great importance. Quantum dot (QD) structures have already been used for 
^ this[8|, |3, llfj making use of the relaxation of two excitons into one bound biexciton on a QD, although QD structures 
cannot achieve the brightness of quantum well (QW) structures. 

The current limitation of QD structures is the low operating temperature of about 4 to 30 K, which is mainly 
due to the decoherence arising from exciton-phonon and hyperfme interactions. QW structures face the additional 
decoherence source due to Coulomb interactions, even at low temperatures. We therefore investigate the possibility 
to use the short time response for the production of entangled photons at time scales for which the biexciton binding 
energy E xx cannot be resolved anymore, i.e. T <C %/E xx , and therefore the Coulomb interactions is negligible. In 
order for this method to be effective, a microcavity is required to extract quickly the entangled photons. In addition, 
1^1 ' the microcavity can be used to tailor the angle dependence of the entanglement. We show that the Pauli exclusion 
principle, which is instantaneous, is sufficient to produce entangled photons from a QW structure with a high yield 
of 10 -2 . This method of producing entangled photons is radically different from the well-known method based on the 



^ . bound biexciton state [8|, |9|, [11|, |lj, [13|, Hi 

The excitation of a QW in a cavity by the external field and emission of the photons due to the radiative recombina- 
£SJ . tion are driven by the interaction of the QW with the photonic modes of the cavity. They are found by quantizing the 
electromagnetic (EM) field in the whole space while taking into account the one-dimensional (ID) spatial modulation 
. of the refractive index n(z) with z axis coinciding with the crystal growth direction. The states of the EM field are 
specified by k — (k||,w, h). Here k|| is the in-plane wavevector, lo is the frequency, h = s,p denotes the polarization 
state. In units H= 1, c = 1 we present the quantized field as (see e.g. Ref. fla) 
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where are the unit polarization vectors (e s lies in the plane of QW), u^(z) is the spatial distribution of the field 
along the z-axis found as the solution of the respective ID scattering problem, p is the coordinate in (x, y) plane 
and aX is the photon creation operator. The summation over k implies the integration over the continuous quantum 
numbers and summation over the discrete ones. 

Entanglement of the photons produced in the course of the radiative relaxation of the pumped semiconductor is 
found considering the two-photon density matrix 

/4i(*)=<»(*)i«i«4«fc«fci»(*», (2) 

where is the state of the semiconductor-photon system. In lowest-order perturbation theory non-vanishing 

terms result from the contribution of the two-photon states into Y$>(t)) , leading to 

4&(*)=^.&(*)*s,a(*) ( 3 ) 
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with g (t) = (0| ag ag where |0) is vacuum of the combined system. The free dynamics in the interaction 

picture is given by 

= (0\a %1 (t)a %2 (t)T + exp j-i J dt'H int (t')\ |0> , (4) 

where a^i) and Hi„ t (t) are the photon annihilation operator and the Hamiltonian of the light- matter interaction, 
respectively, and 7+ is the time-ordering operator. 

The low energy excitations are conveniently accounted by introducing the exciton operators according to \fi) = 
|0), where is the hole-electron pair state either bound or unbound corresponding to energy E^, i.e. Hsc |m) = 
Efj, with Hsc being the Hamiltonian of the nonperturbed semiconductor. For simplicity we assume that the QW 
can be approximated by a 2D plane situated at z = zq. In this case the exciton states are characterized by the spin 
states of the hole and the electron constituting the pair, the center of mass momentum in the plane of the well, K, and 
other quantum numbers, n^, so that = | cr M , s^, K^, n M ). Denoting by ^(x, x') the exciton (hole-electron) wave 
function corresponding to the state \x we represent the exciton operator as = J dx dx! <^(x, x')c S(l {~x!)v a ^ (x), where 
c s and v a are the annihilation operators for the electrons with the spin s and the holes with the spin tr, respectively. In 
terms of the exciton operators the light-matter interaction Hamiltonian is Hi nt = ^2 (A^B^ + A^By) with „4 M and 
An defined as convolutions of A(x) with d CT s <^(x, x) and its conjugate, respectively. Here d nn ' = e/m(n\p\n') 
with (n\ p \n) is the matrix element of the momentum operator between the bands n and n . 

In order to distinguish between the processes of the excitation of the semiconductor and the radiative exciton 
recombination we separate the contributions of the quantized field of emitted photons, A^'(x), and of the classical 
pumping field, A( c/ )(x), into the external field A(x) = A( ? )(x) + A^ c ^(x). In lowest-order perturbation theory we 
can neglect the processes of re-emission and reabsorption, thus obtaining 

H mt *Y,(A q)B »+ A l cl)B l)> ( 5 ) 

which we use in Eq. ([4]). Since in the low-intensity regime the higher-order terms of the light-matter interaction that 
affect the form of the photonic modes in the expression for the density matrix are negligible, the pump field is found 
by solving the respective initial value problem for the cavity alone without the QW. From now on we will omit the 
upper index for the quantized field. 

Expanding the exponential term in Eq. (j4| we obtain various terms depending on ordering of the operators B and 
B* . We leave only the terms describing the response along the directions different from the direction of the incident 
excitation field, where it is not blurred by the non-scattered field and by the linear (single-photon) response. Using 
the 5-functional approximation for the z dependence of the exciton wave function the two-photon amplitude of the 
outgoing photons can be presented as 

%i,fc2^ = u kS z ^ u k2^ e kS z ^ ' ' e S 2 ( z o), (6) 

where the tensor Mk 1 ,k 2 (i) depends only on the direction of propagation of the outgoing photons but neither on their 

polarizations nor on the choice of the solutions of the scattering problems [see the supplement for the explicit form of 
< — > 

■^ki,k 2 (i)J- This information and the effect of the photonic density of states modified by the cavity are contained in 
the amplitudes (zo). If the field distribution inside the cavity has the maximum near Zq this results in accordingly 
amplified two-particle amplitudes ^ . On the contrary, if for one polarization, s or p, the amplitude is significantly 
smaller comparing to the other this can be easily shown to imply significant decrease of the photon entanglement and 
so on. 

The important property of the two-photon amplitudes follows from the conservation of the total in-plane momentum. 
Assuming the normal incidence of the excitation field this leads to the important restriction on the two-photon 
amplitudes ^ cx ^(k^ii +k 2i ||). Since the QW is invariant only with respect to in-plane translations the restrictions 
is imposed only on the in-plane component of the wave vectors of the outgoing photons. 

As Eq. ^ shows, the perturbation theory produces the density matrix corresponding to a pure state of the two 
photon system. Thus we can directly apply the standard machinery for evaluating entanglement of two photons 
as the von Neumann entropy of the reduced density matrix. The dependence of entanglement on the direction of 
propagation of the outgoing photons is very complex. In order to consider the main features of entanglement of 
the emitted photons, we discuss the details of the case of a single QW without the cavity, leaving the details of the 
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case including cavity to another paper. In the case without cavity one obtains a compact expression for the reduced 
single-photon density matrix 

pi(t;k u k 2 ) = e • M kl , k2 (t)(T - e 2 ® e 2 )5?£ ijka (t) • e', (7) 

where 1 is the unit tensor and the argument of p\ shows the dependence of the reduced density matrix on the wave 
vectors of the pair of the outgoing photons. The eigenvalues of the reduced density matrix determine entanglement 
as En = — /5ilog 2 (/5i) — p2log 2 (p2) with pi^ = p\,ij{p\ + p-i)- There are two terms in M^^^i), which describe 
different physical origins of entanglement of the emitted photons. One term describes the creation of entanglement 
due to the Pauli exclusion principle while another term accounts for the effect of the Coulomb interaction. Staying 
in the ultrashort time limit, we can neglect the effect of the Coulomb interaction. 

For the resonance with heavy-hole excitons the non-normalized eigenvalues of the single-particle density matrix are 
proportional to the solutions of 



p 2 - | {cos(2 X ) sin 2 (/3) sin 4 (0) + [l + cos 2 (^)] 



P 
2 

x [l + cos 2 (0)] 2 }+sin 4 (/3)cos 4 (0) = 0, (8) 

where (3 and x are the polar and azimuthal angles on the Poincare sphere describing the polarization state of the 
excitation field, 04 = e ix/2 cos(/3/2) and = e^'*/ 2 sin(/3/2). So that (3 = 0,7r,7r/2 correspond to the left and 
right circular and linear polarizations, respectively, and x/2 is the angle between the axis of the ellipse of polarization 
and the projection of ki on the plane of the QW. Since both excitons have the same energy, the emission has a 
maximum at symmetric angles, i.e. when 9\ = 62 (see Fig. QJi). 

The dependence of entanglement on the direction of the outgoing photons and on the parameters of the excitation 
pulse is relatively simple in this case. Generally, the maximum entanglement is reached near 9 = for the linear 
polarization where the non-normalized eigenvalues are pi.2 — [1 ± cos(/3)] 2 . For linearly polarized excitation field, 
(3 = 7r/2, the pairs are maximally entangled En = 1 while with decreasing the degree of ellipticity entanglement 
decreases. As the special feature the flat maximum near 9 = should be emphasized. There dE/89 cx sin (0) 
for arbitrary polarization of the excitation field. Moreover, for linearly polarized light when x = 7r/2 entanglement 
reaches the maximum, E^ = 1, and is independent of 9 (see Fig. [TJd). It should be noted, however, that as follows 
from Eq. ([5]) the eigenvalues of the reduced density matrix are p\ 2 cx cos 2 (#) in this case, so the signal vanishes in 
the direction 9 = tt/2. 




FIG. 1: (a) The emission directions. The symmetric case corresponds to the heavy-heavy and light-light two-exciton resonances. 
The asymmetric directions 62 = 0± correspond to the heavy-light resonance, (b) - (d) The Dependence of entanglement (vertical 
axes, scale from to 1). (b) En{0,/3 = n/2,x) in the vicinity of the heavy-hole exciton resonance, (c) En{0,/3 = n/2,x = 0) 
in the presence of the cavity with 20 layers of the same optical width at 9 = 0, the index contrast is 1.5, the heavy-hole exciton 
resonance is tuned at the middle of the first stop-band at 9 — 0. (d) En (9, /3,x = 0) near the light-hole exciton resonance. 
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Using Eq. ([6]) in the polarization basis the states are written as \^>) = Ylee' ^ r «>«' l e ' e ')' where 'fee' are found 
convoluting the polarization vectors of the outgoing photons e and e' with ^ . Along the direction 9 m 0, where 
En = 1 can be reached, we obtain 

|*) cx -e* [1 - cos(/3)] |+) |+) - e"'* [1 + cos(/3)] |-) |-) , (9) 

where |+) (|— )) is the state of a photon that is right (left) circularly polarized. As 9 is increased, the entanglement 
is reduced down except for the case (3 = it/2, \ — 7T /2, where the state is of the same structure as for 9 = 0, i.e. 
|<J/) cx |+) |+) — |— ) |— ), as long as 8 < n/2. Along 9 = n/2 the two-photon state is \^f) cx [cos(x)— i cos(/3) sin(x)] |s) |s), 
i.e. the two-photon state is completely disentangled. Here the two-photon state is expressed in terms of the s and 
p polarization cigcnstatcs. Varying the ellipticity of the incoming pulse away from linear polarization (3 = n/2, the 
entanglement is monotonously reduced down to En = for (3 = 0. There the two-photon state reads \^} = |6*i) \9i) 
with \9i) = (i \s) + cos(#) |p)) / \Jl + cos 2 (9), which is also completely disentangled. We would like to emphasize that 
all these states are triplet, that is transform according to the 3D representation of the rotation group. 

Our calculations including the microcavity show that entanglement near 9 = is not affected by the cavity. However, 
in general its angular dependence becomes highly non-trivial, approximately following that of transmittivities, as is 
illustrated in Fig. QJ:. These properties are very useful for tailoring the angle dependence of the emission of the 
entangled photons. The details of the cavity model will be discussed elsewhere. 

For the photons in resonance with the light-hole excitons (lui = uji = E{) the directional dependence of entanglement 
is more complex. The reason is the interaction of obliquely propagating p polarized photons with the light hole 
excitons with the zero projection of the total spin. The characteristic equation in this case differs from Eq. ([5]) by the 
term X {X [p - 1 + sin 2 (/3) sin 2 x] + 2 cos(x) cos 2 (6l) [p - sin 2 (/3)] } in the right-hand side with X = 8 sin(/3) sin 2 (6>). 
This term becomes important for oblique directions thus yielding a richer structure of the angular dependence of 
entanglement, as shown in Fig. [TJi. 

For the photons in resonance with heavy- and light-hole excitons, corresponding to 9\ ^ 82, there is only one 
resonance, when the two-exciton state is made of light-hole and heavy-hole excitons. Respectively, there are only two 
resonant directions with 8\ ^ 82 with 

sin(0 ± ) = sin(0 1 )(^) ±1 . (10) 

The reduced single-photon density matrix is determined by Eq. (O with M ^ cx A + A- ^ 1 — e z <g) e z ^j . Its non- 
normalized eigenvalues are 1 and cos 2 (#). As a result, entanglement monotonously decreases from 1 to as 9 changes 
from to 7r/2 while the direction of the detection of the second photon is determined by Eq. (fT0| . For sufficiently large 
detection angles 9 such that sin(#) > A^/Ei only one resonant direction corresponding to 0_ remains. Interestingly, 
in the asymmetric case entanglement is independent of the the polarization of the excitation field, however, j., 
vanishes if the excitation pulse is circularly polarized. 

Studying entanglement would not be complete without considering yield, which is defined as the ratio of the energy 
flux carried by the entangled pairs of the photons to the flux of the excitation field. For practical purposes it is more 
convenient to use an alternative definition Y = N out /N in , where N out and N in ~ ($/M7) 2 are the number of outgoing 
and incoming pairs, respectively, with $ being the total flux of the external field and £1 being its frequency in the 
stationary frame. Calculating the trace of the single-particle density matrix over the polarization quantum numbers 
we find (in SI units) 

A^E ,(T\Ql_\ 2 



3 7 3 S 

where E x is the exciton energy, pi, 2 are the non-normalized eigenvalues of the density matrix defined by the equations 
studied above, 5* is the area of the excitation spot, T is the duration of the excitation pulse, Q is the common 
interband dipole moment, and eo is vacuum permittivity. The inverse dependence on the pump area is clear since the 
dense excitation more pronounced is the effect of the exclusion principle. The dependence on the pulse duration is 
the consequence of the semiconductor response determined by the polarizations of the linear response rather than by 
the energy input. Substituting the values typical for GaAs, 7 = 1.5 meV, M7 1=3 E x = 1.5 cV, Q = 1.3 • 10~ 24 kg • m/s 



(see e.g. Ref. Il8l ) and using T = 100 fs, S = 7r(20) 2 /1m 2 we find Y ~ 0.02. Such high value of yield is the result of 
the resonant transitions between the many-particle states. 

In conclusion, we have studied the basic mechanism of emission entangled photons by a semiconductor quantum 
well excited by a short pulse. We have developed a "kinematic" theory accounting the effect of the exclusion principle. 
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The dependence of entanglement on the detection angle and on the polarization of the external field is shown to be 
highly nontrivial. We have estimated yield of the considered process and have found it to be rather high owing to the 
resonant transitions between different states. 
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SUPPLEMENTAL INFORMATION 



In the lowest non- vanishing order of the perturbation theory Since the tensor £ (t) does not depend on the 
polarization states of the outgoing photons we omit the hats over fc's and write it solely in terms of the wave vectors 
of the outgoing photons 



t t 

2tt 



M klik2 (t) = ^(0)^(0) dhe-^-"*^-^ / dt 2 e- lE ^~^ 



tH! ^(K w -kx,||)5(K M2 -k^d^,^ 

+e --*x(* I -*=) 5 (K /l2 -k 1 ,||)5(K /Jl -k 2 ,||)d^ 2 ^ 2 Sd^J Q^ 2 (h), 

(12) 

where <g> denotes the tensor product so that in a particular coordinate system (d ® d')ij = did'j, Uk = k is the photon 
energy and are the energies of the hole-electron pair states. The effect of the semiconductor nonlinear response is 
described by the function ^ 1)( i 2 (t), which we present as a sum of the instantaneous term and the memory term 

t 

G^M = {Ba l B^Bl 3 Bl A )P^\t)P^{t) + i J dt> (B^B, 2 e- iH ^-^Di 3tUi )p(l\t')P^\t'), (13) 

o 

where (l9| D\ 3 „ 4 = [Bl , [i?J 4 , H]] and P„{t) are the exciton polarizations of the linear response created by the action 
of the external (classical) field 

t 

pW(t) = -i [ dt'e-^-^A^it'). (14) 



o 
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Neglecting the memory term (or the contribution of the Coulomb interaction) the tensor Mk 1 ,k 2 (^) is directly 
expressed in terms of A + and A_, the amplitudes of the left circular and right circular components of the excitation 
pulse, respectively, 



Mk 1)k2 (t) 



(2^) 5 



5{k Xt \\ +k 2 ,||) (A 2 _e + (g>e + + A 2 + e- ®e_ 



W hh (t) + -Wu(t) 
y 



+ - (e+ ® S_ + e_ ® e+) [W w (i) + W ;/t (t)] 

g 

- -e z ®e z A+A-Wu{t). 



(15) 



Here the indices h and Z denote \a\ — 3/2 and |<r| = 1/2, respectively. The time dependence is described by 
,cr 2 Ki,^ + w ai ^ 2 (u>2, 0J\', t)] with 



II 'n 



(w 1 ,W 2 ',t)= / dtie-*^ 14 ^)^-*^ / <tt 2e -H<>>l+*l.)<.ti-t2) e -iE.mt 



(16) 



where 



E <Jl and 



the exciton-light interaction is quantified by I <71 



|Q| 4 CTl (O)0 ff2 (O) 



_E f7l + _£/ (72 denote single-exciton and two-exciton energies, respectively. The intensity of 

2 

/dq |^ CTl (q)| 2 |^ CT2 (q)| 2 with ^ ff (q) being 

the Fourier transform of the exciton wave function. Here and in Eq. (| 1 5[) we have taken into account the structure 
of the valence band in the semiconductors with the point symmetry and have introduced the common interband 
dipole moment Q = —ie/m (X\p x \Ti) (see e.g. Ref. Il8l). 

The important result immediately following from Eq. (|15|) is that if the external excitation field is circularly polarized 
then the emitted photons are disentangled. Indeed, in this case the tensor M^,^ (t) is represented as a tensor product 
and, hence, so is the single-particle density matrix, i.e. it corresponds to a pure state. 

Representation (|15)) explicitly shows the tensor Mk 1 ,k 2 (t) as a superposition of the amplitudes of the radiative 
decay of different two-exciton states through two channels into the two-photon states. Due to the presence of different 
characteristic frequencies the time dependence of the total amplitude has a complex form, especially at the transitional 
regime with the duration proportional to invcrsed heavy-hole light-hole splitting. However, there are several resonances 
whose amplitudes increases with time and which define the long time response and respectively the time dependence 
of entanglement. In order to extract these resonances we consider the long time limit in the spirit of the Wigner- 
Weisskopf approximation [2fj| extending the limits of integrations over time to infinity. This yields 



w(t) <x e - iE **H [E xx (l + a-b)~ E x (l + a)} 6 + cj 2 - E x 



(17) 



Here we have taken into account the momentum selection rule and have introduced a = sin(#i)/ sin(02) and b stands 
for either 1 or a for w(u>i,u)2) and w(u)2,u>i), respectively, that is depending on particular exciton-photon channel. 
As follows from Eq. (jTTJ) only such terms in Eq. (fT5")) contribute into the long limit which satisfy the special resonant 
condition. The total energy of the emitted pair must be equal to the energy of the two-exciton state. Additionally 
there is the special "kinematic" requirement imposed on the energies of the involved single- and two-exciton states. 



